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ABSTRACT 

We investigate the radio properties of a complete sample of nearby, massive, X-ray bright 
elliptical and SO galaxies. Our sample contains 18 galaxies with ROSAT All-Sky Survey X- 
ray fluxes Fx, 0.1-2. 4kcV > 3 x 10~ 12 ergs _1 cm~ 2 , within a distance of 100 Mpc. For 
these galaxies, we have complete (18/18) VLA radio and Chandra X-ray coverage. Nuclear 
radio emission is detected from 17/18 of the galaxies. Ten of the galaxies exhibit extended 
radio emission; of these ten, all but one also exhibit clear evidence of interaction of the radio 
source with the surrounding, X-ray emitting gas. Among the seven galaxies with unresolved 
radio sources, one has clear, and one has small, cavity-like features in the Chandra X-ray 
images; a third has a disturbed X-ray morphology. Using a radio luminosity limit equivalent 
to £1.4 GHz > 10 23 W Hz -1 to calculate the radio-loud fraction, we find that this misses the 
majority of the radio detected galaxies in the sample. We determine integrated radio-to-X-ray 
flux ratios for the galaxies, QTZx, which are shown to span a large range (factor of 100). We 
calculate the mass-weighted cooling times within 1 kpc, and find hints for an anticorrelation 
with the radio luminosity. We also calculate limits on k/ f, where k is the ratio of the total 
particle energy to that of relativistic electrons radiating in the range 10 MHz-10 GHz and / 
is the volume filling factor of the plasma in the cavity. The k/f distribution is also broad, 
reflecting previous results for larger galaxy clusters. Lowering the X-ray flux limit, at the 
expense of less complete VLA and Chandra coverage, increases the size of our sample to 42 
galaxies. Nuclear radio activity is detected in at least 34/42 of this extended sample. 

Key words: galaxies: active, galaxies: lenticular and elliptical, cD, galaxies: jets 



1 INTRODUCTION 

Super-massive black holes (SMBH) appea r to be common in 
massive galaxies (see e.g.lRichstonej 1 19981 : iFerrarese & Merrittl 
120001 : iKormendvl 1200 ll: iFerraresd I2OO20 . A question of signifi- 
cant interest is the level to which these SMBHs are active i.e. 
their emitted power as a frac t ion of accretion rate or Edding- 
ton luminosity jBest etal]20 05; Nagar et al. 2005; Chiaberg eet al.1 
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20051: i Balmaver de & Capetti |200€ ; ICapetti & Balmaverdd l200fj| : 
Filho et alj|2006l :l Gallo et alj|2008h . Generally power is released 



from active galactic nuclei (AGN) in either radiative or kinetic 
forms. Whereas the radiated power is straightforward to measure, 
using observations across the electromagnetic spectrum, determin- 
ing the kinetic power of AGN is more challenging, particularly for 
low-to-moderate mass galaxies and SMBHs. For the highest mass 
galaxies, however, and in particular those at the centers of clus- 
ters and groups, nature provides powerful, additional tools to probe 
AGN and their impact on their environments. 

Within massive elliptical galaxies, groups and clusters, em- 
bedded AGN have been shown to interact strongly with the sur- 
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rounding hot X-ray emitting medium. This interactio n causes dis- 
turba n ces seen as shocks, ripples, and cavities (e.g. lFabian et al.l 
|2000L l200l 120061: lFormanet al1l2005h . The cavities, which are 
filled with low-density, radio-emitting, relativistic plasma, ap- 
pear as X-ray surface brightness depressions - a consequence of 
the X-ray luminous, thermal gas being displaced by the jets of 
the AGN. X-ray cavities, filled with radio emitting plasma, are 
clearly associate d with the central AGN of nearby clusters such 
as Pe rseus (e.g. 
l2006h . Vir go (e.g. 



Bohringer et al .1 1 19931; iFabian et al.l|2000L [20031 
„. Churazov et al.ll200ll ; lForman et alj|2005f), Cen- 
tauru s (e.g. lTavlor et al. 20061) and Hydra A (e.g. lMcNamara et al.l 
120001) . 

Crucially, studies of cavities allow the total kinetic energy and 
kinetic powe r injected by AGN into their surroundings to be es- 
timat e d (e.g.lChurazov et alj|2002l; lAllen et alj|2 006; Birz an et all 



|2004 iRaffertv et alj|2006l ; iDunn & Fabianll200d l2008h . This m- 
jection is energetically sufficient to prevent catastrophic cooling of 
the X-ray emitting gas, providing a potential solution to the puz- 
zle of why most nearby, massive galax ie s appear 'red and dead ' 
dBower et alj|2006l : ICroton et alj|2006al lbl: ISomerville etai]|2008h . 
For a revi ew of the effect SMB Hs have on their host galaxy's evo- 
lutio n seelCattaneo et al] J2009T). 

Dunn & Fabian (2006, 2008) showed that almost all galaxies 
at the centers of nearby clusters for which the central cooling time 
of the X-ray emitting gas is short (i.e. t coo \ « tu), and for which 
this gas would otherwise cool rapidly and form stars, host power- 
fu l, radio-emitti ng AGN and X-ray cavities, this extended the work 
of lBurnsI d 199Ch who showed that 70 per cent of dominant galaxies 
in rich clusters with cooling cores are rad io loud. Using a larger 
sample of clusters and groups, ISunl d2009l> argue the inverse, that 
all brightest cluster galaxies with radio AGN have cool cores. As 
radio-emitting AGN are seen in almost all clusters with short cool- 
ing times, the duty cycle of these AGN must be high. 

It is important to recognise, however, that the radio sources 
at the centres of clusters are in extreme environments. To assess 
the impact of such processes on the general population of mas- 
sive galaxies, it is essential to extend this work to the general 
populati on of massive e l lipticals. Recent studi es have used the 
' toughton et"al]|2002l) and the NVSS 
exam ine the fraction of galaxies with radio- 
loud AGN. Best et al.l ( 120051) combine these two surveys, and se- 
lect galaxies with redshifts between 0.03 < z < 0.1. They show 
that the fraction of galaxies with a radio-loud AGN depends on 
the stellar mass of the galaxy; the maximum fraction is between 
30 and 40 per cent, significantly below that inferred for nearby 
clusters. (Note that the 'maximum' fraction is obtained for the 
least restrictive constraint on the radio luminosity and the largest 
galaxies in that study.) Other investigations into the fraction of ac- 
tively accreting SMBHs in the nearby Universe have led to sim- 
ilar results, though do not rule out the possibility of significantly 
larger fractions fo r lower radio luminos it y thresholds (see e.g. 
Nagar et alj 120051: Iciiiaberge et alj 120051; I Balmaver de & Capettil 
20061 : ICapetti & Balmaverdel2006l : lFilho et alj|2006h. 



population 01 massive enip 
SDSSH dYork et al.l|2000l ; Isti 
dCondon et alj 199a) to exam 



A nother angle has been taken by Goulding & AlexandeJ 
d2009h who select very local (D < 15 Mpc) galaxies which are 
bolometrically luminous using infra-red observations. In this vol- 



1 Sloan Digital Sky Survey 

2 NRAO (National Radio Astronomy Observatory) VLA (Very Large Ar- 
ray) Sky Survey 



ume limited sample (94 per cent complete) of the 64 galaxies, 17 
host an AGN (27 per cent). 

In this paper we use high quality X-ray and radio observa- 
tions to determine the fraction of active AGN in a complete sample 
of nearby, massive elliptical galaxies, selected to have high optical 
and X-ray fluxes and lie within a distance of lOOMpc. Our sample 
includes both dominant galaxies at the centres of groups and clus- 
ters and field ellipticals with their own halos of hot X-ray gas. We 
show that the fraction of galaxies that are active, i.e. exhibit radio 
emission and/or clear X-ray cavities associated with a central AGN 
is very high: > 90 per cent. We also investigate the prevalence of 
cool cores within this sample. 

Our discussion will proceed as follows. We begin in Section 
[2]with a description of the sample and selection criteria. In Section 
[3] we describe the radio data selection and processing. In Section 
[4] we describe the X-ray data reduction and preparation. In Section 
[5] we discuss individual sources and in Section [6] we conduct a 
population study. Section|7]extends some aspects of the analysis to 
a larger sample of galaxies, at the expense of incomplete VLA and 
Chandra coverage. We conclude in Section|9] 

Throughout this paper, we assume a standard flat ACDM cos- 
mology with _ffo=70 km s _1 Mpc -1 and f2 m = 0.3. 



2 SAMPLE SELECTION 

Our study is focused on a clearly defined and essentially complete 
sample of nearby elliptical galaxies with measured o ptical, X-ray 
and ra dio properties. We start from the catalogue of iBeuing et alj 
( 1999). The northern part o f this catalogue con tains all E and E/S0 
galaxies from the study of iFaber et all dl989h with 5^0° and 
magnit udes brighter than Bt = 13.5 mag. T he southern pa r t of th e 
Beuin sTet alj dl999h catalogue is drawn from lBender et alj dl989h : 
it contains all E and E/SO ga laxies with S < 0° and Bt < 13.5 
mag from (Fab er etal] 1 19891), as well as galaxies from the ESO 
Lauberts-Valentijn Catalogue dLauberts & Valentiinl 19891) brighter 
than Bt = 13.5 mag and with morphological type T id or T ncw ^ 
—3. In addition, the catalogue includes galaxies with 5^0°, 
Bt < 13.7 and numerical class T ^ —3 from the Third Ref- 
erence Catalogue for Bright Galaxies dde Vaucouleurs et alJI 199lb 
and the Tully Nearb y Galaxies Ca t alogue dTullvll988H 3 H 4 l 

The sample of Beui ng et alj d 19991) contains 530 early type 
galaxies; 313 in the southern and 217 in the northern parts. The 
distribution of galaxies of different types is as follows: 49 per 
cent ellipticals (T = —5), one of which is a compact elliptical 
(T = —6); 15 per cent cD-type galaxies (T = —4); 26 per cent 
E/SOa (T = -3); 7 per cent SOs (T = -2); and 3 per cent 
later types dde Vaucouleurs et alj 1991 ). From their magnitude- 
distribution function, Beuin g et al l dl999h judge that their sample 
starts to become incomplete at Bt ~ 13.0. However they state that 
it is still 90 per cent complete at Bt = 13.5, assuming a spatially 
hom ogeneous distributio n. 

IBeuing etal.ldl999h provide X-ray luminosities, or upper lim- 
its on the X-ray luminosity, for 293 galaxies in their sample, based 



3 Beuing etal. (1999) adopt a 


de Vaucouleurs et al. 


d 1 99 lh and 


de Vaucouleurs et al. 


Jl99ll) found 



fainte r Bt threshold for the 
iTuiivl i 19881) catalogues, since 
hints of a zero-point offset of 
0.12 - 0.18 mag. 

4 The numerical class constraint excludes at least one well known X-ray 
bright galaxy, NGC 1275, as it is classed as a peculiar galaxy, T = 99 
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on data from the ROSAlffl All-Sky Survey (RASS). They exclude 
galaxies which are not at the centre of the bright X-ray emission, 
but embedded in it, unless they clearly stand out against it. The up- 
per limits obtained for some galaxies were classed as meaningless 
if they resulted from short exposures, or if the galaxy has a distance 
> Mpc. Their final sample consists of all normal galaxies 

showing individual emission, as well as galaxies at the centres of 
clusters and groups, if the X-ray emis sion is centred on the galaxy 
and reasonably symmetric around it. |0' Sullivan et alj d200ll) up- 
dated these Z/x values, also incorporating pointed, follow-up X-ray 
observations made with the ROSAT Position Sensitive Proportional 
Counter (PSPC), where a vailable. For our s ample selection, we use 
the X-ray luminosities of O'Sull ivan et alj 1 2001 ), converted to our 
reference cosmology. 207 galaxies in the lO' Sullivan et "all feOOlh 
study have X-ray luminosity measurements as opposed to upper 
limits. 

From the initial catalogue of optically bright galaxies with X- 
ray luminosity measurements, we identify all galaxies with 0.1 — 
2.4keV fluxes > 3 x 10~ 12 erg s _1 cm -2 . We also impose an 
upper distance limit of 100 Mpc. Together, these two criteria en- 
sure that all galaxies studied should be sufficiently X-ray bright 
and nearby to allow the properties of their diffuse, X-ray emitting 
gas to be studied in detail with Chandra. 

We note t hat the X-ray luminosities reported by 
lO' Sullivan etaT] ( 1200 lh are total luminosities, including the 
contributions from X-ray binaries, AGN and other compact 
sources, as well as the hot, gaseous components. For the largest 
elliptical galaxies the contributions from b i naries is expected to be 
minimal (see e.g. Kim & Fabbianol 120041 : [Humphrey et alj 1200 8t 
iBrough et al]|2008ir As we show in this paper, the impact of the 
central AGN on the total X-ray emission is also typically low. Only 
in ~ 10 per cent (2/18) of our galaxies does the AGN outshine the 



diffuse gas at X-ray wavelengths (see Sections |5land |6. 1| 6 | For the 
other galaxies, the X-ray emission from the central AGN typically 
accounts for less than 1 per cent of the total. 

The X-ray properties of the parent sample and X-ray flux and 
distance cuts are shown in Fig. [JJ The resulting sample consists 
of 18 galaxies, which are detailed in Table [TJ VLA radio data and 
Chandra X-ray observations are available for all 18 targets. 

Two of the eighteen sources, IC310 and NGC4203, are 
anomalous in that their X-ray emission is dominated by the central 
agnO These two galaxies are the two most late type (SO) galax- 
ies in the sample (Table [TJ. For these reasons, although IC 310 and 
NGC 4203 formally meet the optical, X-ray and distance selection 
criteria, we flag or exclude them from certain population studies. 



3 RADIO DATA 

The radio data for the targets were obtained from the NRAO VLA 
archivefj Multiple observations of each source are available, for 



5 Rontgen Satellit (ROSAT) 

6 Two of the eighteen sources, IC310 and NGC 4203, are anomalous in 
that their X-ray emission is dominated by their central AGN. These two 
galaxies are also the two most late-type (SO) galaxies in the sample (see 
TablefJJ. For these reasons, although IC 310 and NGC 4203 formally meet 
the optical, X-ray and distance selection criteria, we flag or exclude them 
from certain population studies. 

7 The X-ray luminosities of O'Sullivan et al. 12001). include the contribu- 
tions of point sources in the galaxies. 

8 The VLA (Very Large Array) is operated by the National Radio Astron- 
omy Observatory (NRAO). The National Radio Astronomy Observatory 
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Figure 1. The dis tance-luminosity distribution of the galaxy sample of 
iBeuine et alj . 1999 with the flux and distance cuts used to obtain our final 
sample. The 1 8 galaxies in our sample are highlighted in red. Galaxies for 
which only upper limits on the X-ray luminosity are available are shown 
with arrows. IC310 and NGC 4203, for which the X-ray emission is dom- 
inated by the central AGN, are highlighte d with blue ci rcles. The X-ray 
luminosities shown are taken from O'Sullivan et al. 1 200 J) but converted to 
H = TOkms" 1 Mpc" 1 . 



various wavelengths and configurations. For each source we select 
and present the radio map best suited to the present purposes. 

In selecting the data to be analysed, preference was given 
to more recent observations, and to those with time on source 
^ 5 minutes. As steep spectrum emission is brighter at lower fre- 
quencies, preference was also given to observations performed at 
1.5 GHz. In addition, A-configuration observations at 1.5 GHz, or 
B-configuration observations at 5 GHz were desirable in order to 
provide arcsecond resolution for comparison with X-ray images 
from Chandra. 

The data wer e reduced in the standard manner using 
After an initial editing of the data, absolute 
amplitude and phase calibration were performed on each dataset 
using the scripts VLAPROCS and VLARUN. For datasets in which 
the flux calibrator was resolved, a model was used (if available) for 
the calibration. If bad data were still present after the initial cali- 
bration, those data were flagged and the calibration was repeated. 

An image of each source was created using the task IMAGR. 
Side-lobes from outlying sources were removed by using multi- 
ple facets while imaging. Proper placement of the facets was de- 
termined using the task SETFC, which was set up to search a 0.5 
degree radius for sources in the NVSS catalogue with flux ^ 10 
mjy. Typically, sources outside of that range are either not bright 
enough or too far from the pointing centre to have an appreciable 
effect on the quality of the image. In any event in which SETFC 
failed to include significant sources, facet positions were set manu- 



is operated by Associated Universities, Inc., under cooperative agreement 
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Table 1. The Sample 



Source 


Alternate Name 


RA a 


Dec a 


Redshift 3 


Distance 


Type 6 


log(L x ) 






J2000 


J2000 




Mpc 




erg s 



IC310 


— 


03h 16m 43.00s 


+4 Id 19m 29.4s 


0.018940 


67.92 


-2.0 


- 


42.60 


IC 1860 


— 


02h 49m 33.7s 


-31dllm21s 


0.022902 


96.59 


-4.7 


BCG 


42.77 


NGC499 


— 


Olh 23m 11.5s 


+33d27m38s 


0.014673 


59.15 


-2.8 


- 


42.35 


NGC507 


— 


Olh 23m 40.0s 


+33dl5m20s 


0.016458 


71.99 


-3.2 


- 


42.76 


NGC533 


— 


01h25m31.36s 


+01d45m32.8s 


0.018509 


68.23 


-4.8 


- 


42.29 


NGC708 


A262 


Olh 52m 46.48s 


+36d 09m 06.6s 


0.016195 


59.15 


-4.8 


BCG 


43.09 


NGC 1399 




03h 38m 29.08s 


-35d 27m 02.7s 


0.004753 


19.40 


-4.2 


BCG 


41.69 


NGC 1404 




03h 38m 5 1.92s 


-35d 35m 39.8s 


0.006494 


19.40 


-4.7 




41.25 


NGC 1550 




04h 19m 37.93s 


+02d 24m 35.7s 


0.012389 


51.95 


-3.9 




42.86 


NGC 4203 




12h 15m 05.06s 


+33d 11m 50.38s 


0.003623 


17.38 


-2.7 




41.24 


NGC 4406 


M86 


12h 26m 11.74s 


+12d 56m 46.40s 


-0.000814 


17.06 


-4.7 




42.11 


NGC 4472 


M49 


12h 29m 46.76s 


+08d00m01.71s 


0.003326 


17.06 


-4.7 




41.49 


NGC 4486 


M87 


12h 30m 49.4s 


+12d 23m 28s 


0.004360 


17.06 


-4.3 


BCG 


43.01 


NGC 4636 




12h42m49.9s 


+02d41ml6s 


0.003129 


17.06 


-4.8 




41.65 


NGC 4649 


M60 


12h 43m 39.66s 


+1 Id 33m 09.4s 


0.003726 


17.06 


-4.6 




41.34 


NGC 4696 


Centaurus Cluster 


12h 48m 49.3s 


-41dl8m40s 


0.009867 


39.65 


-3.9 


BCG 


43.29 


NGC 5044 




13h 15m 23.97s 


-16d23m07.9s 


0.009020 


32.36 


-4.8 




42.80 


NGC 5 846 




15h 06m 29.29s 


+()ld 36m 20.24s 


0.005717 


24.55 


-4.7 




41.71 



The sample of ga l axies. a Positions and redshifts obtained from NASA/IPAC Extragalactic Database (NED). b The Type shows the T-type 
Ide Vaucouleurs et all ll99ll) with E=-5, E/S0=-3, S0=-2, S0a=0. Also shown is whether the galaxy is a Brig htest Cluster Galaxy, (BCG), defined as the 
galaxy (likely) to be the dominant one in a cluster listed in the Abell catalogues iAbelj|l958tlAbell et alj|l989l) . Also listed are the total X-ray luminosities 
(Section^. 



ally. For those sources with sufficient signal-to-noise ratios, imag- 
ing and phase-only self-calibration were then performed iteratively, 
until the theoretical noise was reached or until the quality of the 
map ceased to benefit from the iterations. Observational parame- 
ters and data for the VLA observations are tabulated in Table|2] 

For NGC 1404 the VLA radio map from 2007 shows a faint 
point source coincident with the galaxy centre. However the de- 
tection is only just above the level of the noise in the data, and 
so cannot be determined accurately. The NVSS catalogue also 
lists a source at the location of NGC 1404. We list this detection 
separately in Table [2] During the analysis of the radio data for 
NGC 4472 we found that the data from the A configuration has a 
(it, v) coverage such that the large scale emission was not detected. 
As a result of this, we checked all other A-configuration observa- 
tions for undetected large scale emission using C-configuration data 
present in the VLA Archive. This problem was not found in any 
other galaxy apart from NGC 4472. 

The reduction and an alysis of the 1.5 G Hz radio data for 
NGC 4486 are described bv lHines et ail 1 19891) . 



4 X-RAY DATA REDUCTION 

The Chandra X-ray data were cleaned and reprocessed using the 
CIAO software and latest calibration files (CIAO v4.1.2, CALDB 
v4.1.1). Hot pixels and cosmic ray afterglows were identified and 
excluded. Charge-Transfer Inefficiency was corrected for and stan- 
dard grade selection applied. Light curves were generated from rel- 
atively source-free regions of the detectors and used to screen the 
data for background flares. For observations made in Very Faint 
(VF) mode, the additional information available for identifying and 
screening background events was utilised. Background files were 
generated from the CALDB blank-field data-sets. These were pro- 
cessed in a similar manner to the target observations. The back- 
ground data-sets were normalised by the ratio of the 9 — 12 keV 



count rates in the source and background data-sets. Point-sources 
were identified using the WAVDETECT wavelet-transform proce- 
dure, and excluded. 

To model other, possible background components not ac- 
counted for in the blank-sky fields, background spectra were also 
extracted from relatively source-free regions of the detectors (e.g. 
the ACIS-S 1 chip for ACIS-S data) and compared with blank-sky 
fields. Any excess soft emission components so detected were fit- 
ted with a single temperature MEKAL model with a fixed solar 
abundance and scaled appropriately in the subsequent analysis. 
NGC 5846 exhibits clear excess soft background emission. Note, 
however, that care is required here since in many cases the emis- 
sion from the cluster/group extends onto the S 1 chip. 

For each galaxy, concentric annular regions were selected to 
give approximately constant signal-to-noise in the 0.6 — 7.0 keV 
band. Annuli were centered on the peak of the X-ray emission. 
Where the Chandra data revealed the presence of a central point 
source, the minimum radius was adapted to exclude it. The signal- 
to-noise ratio for the annuli was adjusted iteratively to obtain a 
number of regions ranging from 10 to 25. 

The 0.6 — 7.0 keV spectrum for each spherical shell was 
modelled as a sin gle temperature plasma, using the MEKAL (e.g. 
iMewe et"al]| 19951) model, with a PHABS absorption model fixed to 
the galactic Nn- The spectra for all annuli were modelled simulta- 
neously to determine the deprojected gas temperature, metallicity 
and density profi les for the galaxies. The modified C-statistic avail- 
able in XSPEC ( Arnaud 1996, 2004) was used in all spectral fitting. 



5 NOTES ON INDIVIDUAL SOURCES 

Our sample is essentially unbiased with respect to the morphol- 
ogy of the X-ray emission from the galaxies. X-ray emission from 
the hot, diffuse gaseous halos and any AGN that may be present 
will both contribute to the measured X-ray flux. In this Section, 
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Table 2. VLA Data and Observational Parameters 



Source 


Frequency 


Config 


Date Time on Source 


Peak Flux Total Flux Density 


RMS 


Beam 


PA 




(GHz) 






(sec) 


(Jy/beam) 


(Jy) 


(/jjy/beam) 


(arcsec) 


(deg) 


IC310 


4.885 


C 1984 May 01 


1170 


1.36X10- 1 


1.98X10- 1 


73.9 


3.99x3.93 


-43.55 


IC 1860 


1.365 


A 2007 Jun 08 


8070 


1.23xl0~ 2 


1.39xl0~ 2 


60.2 


2.92 x 1.42 


11.83 


NGC 499 


1.365 


A 2007 Jun 08 


6070 


<1.18xl0- 4 


- 


35.8 


1.43x1.16 


-62.60 


NGC 507 


1.525 


C 1984 May 14 


1770 


8.61 xl0~ 3 


8.45 xl0~ 2 


87.1 


17. 58x 12.11 


-81.43 


NGC 533 


4.885 


C 1997 Jul 05 


230 


5.06xl0~ 3 


8.91X10" 3 


93.3 


5.08x4.60 


-14.47 


NGC 708 


1.365 


A 2000 Dec 13 


1520 


6.57X10" 3 


8.25 xl0~ 2 


55.5 


1 .46 x 1 . 1 8 


-49.36 


NGC 1399 


1.465 


H 1 1983 Dec 16 


27910 


1.89X10- 2 


4.63 xlO" 1 


72.1 


4.15x2.80 


38.46 


NGC 1404 


1.365 


A 2007 Jun 08 


OUOU 


^ v in~4 

j. jU X 1U 




85.8 


3. 64 x 1.48 


10.82 


NGC 1550 


1.365 


C 1997 Sep 18 


410 


2.68xl0~ 3 


1.27xl0~ 2 


76.7 


18.6x16.8 


-84.14 


NGC 4203 


4.885 


A 1999 Sep 05 


423 


1.66X10- 2 


1.80X10- 2 


93.6 


0.48x0.39 


63.12 


NGC 4406 


4.885 


C 1984 Jun 03 


510 


5.09xl0 -4 


4.91 xl0~ 4 


67.1 


4.32x3.86 


55.55 


NGC 4472 


1.489 


C 1986 Nov 26 


192 


1.19X10- 1 


2.23 xlO" 1 


250 


17.19x14.24 


-37.6 


NGC 4486 


1.466 


A 1989 Mar 01 


10020 


1.74x10° 


1.38x10 1 


32.5 


0.80x0.80 


0.00 


NGC 4636 


1.425 


A 2002 Mar 12 


12440 


4.80xl0" 3 


6.45 xlO" 2 


45.8 


1.42x1.33 


-18.42 


NGC 4649 


1.465 


B 1984 Jan 24 


4410 


1.72X10" 2 


2.82xl0~ 2 


29.7 


4.51x3.64 


44.76 


NGC 4696 


1.565 


A 1998 Apr 23 


3635 


3.58X10- 1 


3.64x10 


118 


4.40x2.09 


0.50 


NGC 5044 


1.465 


A 1992 Nov 27 


8560 


2.71xl0~ 2 


2.99 xl0~ 2 


32.8 


1.77x1.03 


-16.83 


NGC 5846 


1.465 


A 2002 Mar 12 


9220 


9.93xl0~ 3 


1.02xl0~ 2 


20.8 


1.32x1.18 


-17.86 


NGC 1404+ 


1.400 


D 1993 Oct 07 




3.42X10" 4 


3.90X10" 4 


450 


A ^ \y A ^ 


U 


1 Hybrid configuration (A and B). t From NVSS,as the point t 


iource in the A-configuration image is only just above the noise. 






Table 3. X-ray data parameters 






















Source 


Date 


ObsID 


Detector Mode 


Exposure 1 


IV H 










IC310 


23 Dec 2004 


5597 


ACIS-I FAINT 


25.2 


1.30 










IC 1860 


12 Sep 2009 


10537 


ACIS-S VFAINT 


37.2 


2.40 










NGC 499 


4 Feb 2009 


10865 


ACIS-S VFAINT 


5.1 


5.30 












5 Feb 2009 


10866 


ACIS-S VFAINT 


8.0 


5.30 












7 Feb 2009 


10867 


ACIS-S VFAINT 


7.0 


5.30 












12 Feb 2009 


10536 


ACIS-S VFAINT 


18.4 


5.30 










NGC 507 


1 1 Oct 2000 


2882 


ACIS-I VFAINT 


43.6 


5.24 










NGC 533 


28 Jul 2002 


2880 


ACIS-S VFAINT 


34.8 


3.10 










NGC 708 


03 Aug 2001 


7921 


ACIS-S VFAINT 


110.6 


5.37 










NGC 1399 


18 Jan 2000 


319 


ACIS-S FAINT 


55.9 


1.34 










NGC 1404 


13 Feb 2003 


2492 


ACIS-S FAINT 


29.0 


1.45 










NGC 1550 


08 Jan 2002 


5800 


ACIS-S VFAINT 


44.5 


11.50 










NGC 4203 


10 Mar 2009 


10535 


ACIS-S VFAINT 


41.4 


1.10 










NGC 4406 


07 Apr 2000 


318 


ACIS-S FAINT 


12.6 


2.62 










NGC 4472 


12 Jun 2000 


321 


ACIS-S VFAINT 


29.6 


1.66 










NGC 4486 


05 Jul 2002 


2707 


ACIS-S FAINT 


98.4 


2.54 










NGC 4636 


26 Jan 2000 


323 


ACIS-S FAINT 


44.2 


1.81 










NGC 4649 


30 Jan 2007 


8182 


ACIS-S VFAINT 


49.2 


2.20 












20 Apr 2000 


785 


ACIS-S VFAINT 


19.4 


2.20 










NGC4696 2 


01 Apr 2004 


4954 


ACIS-S FAINT 


88.0 


8.06 










NGC 5044 


19 Mar 2000 


3225 


ACIS-S VFAINT 


83.1 


4.93 










NGC 5846 


24 May 2000 


788 


ACIS-S FAINT 


22.1 


4.26 







1 The final exposure in ks after all screening. 2 Also analysed in lTavlor et alj ^20061) . t The Galactic column density, JVjj, is in units of 10 20 atomcra 2 
determined from Hi studies IDickev & Loc kman 1990). 



we discuss the X-ray and radio morphologies for the galaxies in 5.1 IC 310 
the sample, comment on the presence of nuclear X-ray emission, 
and examine the evidence for interactions between the central ra- 
dio sources and surrounding X-ray gas. Several of the galaxies in IC310 i s an SO galaxy in the south-west region of the Perseus 
the sample have previously been studied by other authors. cluster. The radio morpho logy is that of a head-tail galaxy (e.g. 

Sijbring & de Bruvnl [T998) with a bright core and extended tail 
of length around 400 kpc. Our data show the radio tail discon- 
nected from the much brighter head. There is only an observation 
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Figure 2. Adaptively smoothed X-ray images (black/purple for faint emission, red/white for bright emission) with radio contours overlaid. From top left to 
bottom right: IC 310, IC 1860, NGC 499, NGC 507, NGC 533, NGC 708, NGC 1399, NGC 1404 and NGC 1550. The Chandra image of NGC507 also shows 
artefacts associated with chip gaps in the ACIS-I array. The radio contours for NGC 1404 are those from the NVSS survey. 



at 4.9 GHz which has sufficient resolution to show the structure of 
the source. 

The X-ray emission is dominated by a central point source, 
coincident with the radio core. The extended X-ray halo is faint, 
reflecting the relatively late-type nature of the galaxy. A small X- 
ray extension is observed in the same direction as the radio tail. The 
current data do not allow any detailed study of interactions between 
the radio and X-ray emitting plasmas. 

5.2 IC 1860 

IC 1860 is the dominant galaxy of the IC 1860 group, which is itself 
pa rt of Abell S3 1. The dominant galaxy of Abell S301 according 
to lHudson et all l l200ll) is the Seyfert 2 galaxy IC 1859, which is 
6.5 arcminutes distant from IC 1860. However, XMM-Newton data 
(Observation 0146510401) shown in Fig.|2] IC 1860 sits at the peak 
of the diffuse X-ray emission, and the position of IC 1859 is off to 
the west. 



The radio data of IC 1860 show a clear point source with faint 
extensions in the northeast and southwest directions. The X-ray 
data show an extension towards the southeast. However there is no 
clear correlation between the faint radio extensions and the X-ray 
gas. 

5.3 NGC 499 

NGC 499 lies near NGC 507 in the NGC 507 galaxy Group. There 
is no radio detection for this galaxy down to a level of 1.18 x 10~ 4 
ly, identifying NGC 499 as abnormally radio quiet for its X-ray 
luminosity. The RMS noise in the VLA data is 35.8 /ily/beam, 
which is fairly typical of the observations presented here (see Table 
O. There appears to be a depression in the X-ray emission to the 
south-east, visible in Fig. [2] However, whet her this is the re s ult of a 
past outburst of AGN activity is not clear. I Kim & Fabbianol dl995h : 
Paolil hTet al.l d2003b report evidence for tidal interactions between 
NGC 499 and NGC 507. 
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Figure 2. (cont.) From top left to bottom right: NGC 4203, NGC 4406, NGC 4472, NGC4486, NGC 4636, NGC 4649, NGC 4969, NGC 5044 and NGC 5846. 



5.4 NGC 507 

NGC 507 is the brightest member of a group of galaxies that is part 
of the Pisces Cluster. Our radio map has comparatively poor spa- 
tial resolution but show that the source is highly extended with two 
bright lobes in a roughly east- west configuration. ROSAT X-ray ob- 
servations showed NGC 507 to have an extended, bright X-ray halo 
with a luminosity comparable to that of poor clusters jKriss et all 
1 19831 : IK m & F abbian oi l 1995b . The Chandra data show some dis- 
turbance at the centre of the cluster with indicatio ns of a shell of en - 
hanced emission around the eastern radio lobe (Kra ft et al. I l2004h . 
However, any signs of interaction are less clear around the western 
lobe, though the secondary western, X- ray bright lump ma y have 
been uplifted by the action of the AGN dForman et alj|200lh . 

5.5 NGC 533 

NGC 533 is the dominant galaxy of a group of the same name. Our 
4.9 GHz VLA radio data show that the radio emission is domi- 
nated by a central point source, although some extension to the 



west and southwest is detected at modest significance. The X-ray 
peak is coincident wit h the central radio source (see S [2009; 
iGastaldello etai]|2007l; IPiffaretti et aTll2005l) . On smaller scales, 
two crescent shaped depressions in the X-ray emission are ob- 
served, either side of the peak, with a roughly northeast-southwest 
orientation (see also Fig.|4j. The resolution of the current radio data 
is insufficient to determine unambiguously whether these features 
are the result of AGN interaction. The cavities are therefore classi- 
fied as 'possible' rather than 'definite'. The larger scale X-ray emis- 
sion also exhibits northeast-southwest elongation, with a bay-like 
feature to the southwest. Low-frequency radio observations would 
be required to investigate whether this bay is associated with a past 
generation of AGN activity. 

5.6 NGC 708 

NGC 708 is the brightest galaxy of the Abell 262 cluster and hosts 
the radio source B2 0149+35. The VLA data reveal a central point 
source and two lobes in an east-west c onfiguration. The X -ray 
emission from the cluster is described bv lBlanton et al.l J2004I) . A 
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The X-ray emission is peaked on the central AGN. Regions of en- 
hanced X-ray emiss ion are observed just beyond the radio lobes 
(Shur kin et all2008h . 




Figure 3. The adaptively smoothed XMM-Newton X-ray emission sur- 
rounding IC 1860 using the same colour scheme as Fig. [2] showing the lo- 
cation of IC 1 859 to the west. The radio contours from the observation of 
IC 1 860 are overlaid. The dominant galaxy, at least according to the X-ray 
emission, appears to be IC 1 860. 




5.8 NGC 1404 

NGC 1404 is a member of the Forna x Cluster, in the proce ss of 
falling towards the centre (NGC 1399, lMachacek et al.|[200^) . Our 
VLA A-configuration data show a faint central point source. The 
radio source is also clearly detected in the NVSS catalogue. The 
X-ray emission map is fairly symmetric and also does not show 
any features that would indicate significa nt interaction between the 
central AGN and surrounding X-ray gas dMachacek et al . 2005). 



5.9 NGC 1550 

NGC 1550, also known as NGC 1551, is the dominant galaxy of the 
NGC 1550 group. The source exhibits an unusual radio morphol- 
ogy with two peaks: one approximately centered on NGC 1550 and 
the other offset ~ 45 arcsec to the west. The group is unusually 
X-ray luminous for its temperature and velocity dispersion. The 
X-ray emission is elliptical, with the semi major axis in th e same 
direction as the axis of the radio source feunet all 12003). A fil- 
ament of X-ray emission traces the eastern radio lobe, wrapping 
around the southern part. The western radio peak appears uncorre- 
lated with the X-ray emission. It is possible that the western radio 
peak is actually due to a background source, although the NASaQ 
Extragalactic Databas^] lists no source at this position. It is also 
possible that radio emission from NGC 1550 has expanded in the 
surrounding ICM in a highly asymmetric way to create the western 
peak. 



5.10 NGC 4203 

NGC 4203 is an SO galaxy hosting a low-luminosity low- 
ionisation nuclear em ission region (LINER) type AGN (LLAGN) 
jlvomotoet"aT] ri998). Point like radio emission is detected in our 
4.9 GHz VLA data. The Chandra observation shows an X-ray 
structure similar to IC310. There is very little diffuse emission, 
with most of the X-rays coming from a point source coincident 
with the radio emission. 



Figure 4. Adaptively smoothed X-ray image of NGC 533 using the same 
colour scheme as Fig. [5] with the radio contours overlaid in black, and 
the two depressions in the X-ray emission mentioned in the text high- 
lighted. Their radii/semi-major axes are 0.88 for the north-eastern and 
0.92 X 1.31 kpc for south-western cavities. 



co mbined study of low frequency radio and X-ray data is presented 
bv lClarkeetai] 12009), showing the interaction of the radio source 
with its surroundings. 



5.7 NGC 1399 

NGC 1399 is the dominant galaxy in the Fornax Cluster (Abell 
S0373). The radio and X-ray properties are described by 
IShurkinetal1 d2008). The radio source exhibits a nucleus with two 
tightly collimated jets that feed into more extended radio lobes. 



5.11 NGC 4406 

Also known as M 86, NGC 4406 is the dominant member of a small 
group falling at high velocity towards the centre of the Virgo clus- 
ter (NGC4486/M 87). Our 4.9 GHz VLA radio data reveal a cen- 
tral point source. XM M-Newton and Chandra o bser vations of the 
galaxy are discussed by Finog uenov et al. (2004) and Rand all et al.l 
( 2008), respectively. Our Chandra image shows only the core of the 
X-ray emitting gas. There are no clear X-ray features indicating an 
interaction between the radio and X-ray plasmas. 



10 National Aeronautics and Space Administration 

11 NASA/IPAC Extragalactic Database (NED) is operated by the Jet 
Propulsion Laboratory, California Institute of Technology, under contract 
with the National Aeronautics and Space Administration. 
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5.12 NGC4472 

Also known as M 49, this Seyfert 2 galaxy is the optically brightest 
elliptical in the Virgo cluster. It is also the dominant member of a 
small group. The cent ral radio source has a typical FRI structure 
but is unusually weak (Ekers & Kotanvill 19781) . We present the data 
from the VLA in C configuration. As noted in Section [3] the data 
from the A configuration has a (u, v) coverage such that the large 
scale emission is not detected, and hence we present the C configu- 
ration data. The VLA data is a comparatively short observation, but 
the morphology matches that seen i n the FIRSTS surve y observa- 
tion, presented in lBiller et all d2004l) ; lAiien et al.U2006l) . There are 
cavities in the X-r ay emitting gas corresponding to the radio lobes 
jBiller et alj|2004> . 



5.13 NGC4486 

NGC 4486, also known as M 87 (Virgo A), is the dominant galaxy 
at the centre of the Virgo Cluster. It hosts the well studied radio- 
bright AGN 3C 274. A bright radio jet is also clearly visible in 
the VLA images, although the contours plotted have been cho- 
sen to highlight the diffuse emission surrounding the core. The jet 
is also clearly seen at X-ray wavelengths. The large scale struc- 
ture of t he ICM and radio emission in NGC 4486 is complex 
(see e.g. iHines et alj|l989l; lOwen et alJl200Ct lYoung et alj | 2002 



Pi Matteo et al.l 120031 ; iForman et alj 120051 : iKovalev et alj 12007 : 
Forman et al]2007h . 



5.14 NGC 4636 

NGC 4636 (also known as NGC 4624) lies in the southern regions 
of the Virgo cluster, some 2.6 Mpc from NGC 4486. The radio data 
show a central core and two clear lobes. The X-ray morphology is 
highly disturbed and suggests that a much larger AGN outburst has 
occurr e d in the past (see also I Jones et alj d2002h : lo' Sullivan et alj 
d2005h : lBaldietal.U2009l) ). 



5.15 NGC 4649 

NGC 4649, also known as M60 is located in a group at the 
eastern edge of the Virgo Cluster. Recent detailed analyses of 
the radio and X-ray p r operties of th i s sour c e can be found 
Randall et al] ( 120041) ; IShurkin et all d2008l) ; iHumphrev etal] 



(2008). In the deeper Chandra data presented here, faint fingers of 
X-ray emission are observed to surround the radio emission, indi- 
cating a channel carved by the AGN jet. Cavities are seen in the 
X-ray gas at positions correlat ed with the radio em ission from the 
AGN, confirming the results of lShurkin et aljj2008h . These appear 
as "wedges" pointing to the surface brightness peak in Fig. [2] To 
clarify the location of the cavities we subtract a spherical profile 
(beta model) which is fitted to the adaptively smoothed image. The 
subtracted image is shown in Fig. [5] 



5.16 NGC 4696 

NGC 4696 is the dominant galaxy of the Centaurus Cluster (Abell 
3526) and hosts the radio source PKS 1245-41. The radio emission 
has two clear lobes. High frequency radio emission shows that the 




Figure 5. The central regions of NGC 4649 where a spherical (beta) model 
subtracted from the adaptively smoothed image. The colour scheme is the 
same as in Fig. [2] The depressions associated with the radio emission (black 
contours) are clearly visible to the north and south of the core. 



jets in itially travel east- west before both turning south dTavlor et alj 
2002). The interaction of the radio source with the surrounding 
intra-cluster medium is s t rikingly clear and is discussed in detail b y 
ISanders & Fabianl ( 120021 ) ; lFabian et"aHd2005l) : lTavlor et alj hoOrj) . 



5.17 NGC 5044 



NGC 5044 is the central member of a rich group dRickes et alj 
l2004h ). Our VLA A-configuration data detect a central point radio 
source. The X-ray e mi ssion at the centre of thi s group is disturbed 
dBuote et alJkfXB al ibi). iGastaldello etal] d2009l) cite evidence for a 
pair of bubbles, one north and one south of the X-ray core. Al- 
though there is no high frequency radio emission associated with 
them, Ha and [ N II ] emissio n does appear to be cor r elated with 
depre ssions in the X-ray gas (Gastaldello et al. 2009; Caon et al. 
2000). More recent GMRlQ observation show that the southern 
cavity is filled by radio plas ma emitting at 610 and 235 MHz radio 
emission dDavid et alj|2009l) . 



5.18 NGC 5846 

NGC 5846 is one of a dominant galaxy pair (with NGC 5850) of the 
NGC 5846 group. Our VLA A-configuration data detect a central 
point radio source. The morphology of the X-ray emission is highly 
suggestive of an AGN interaction, with a clear circula r cavity to the 
north-east of the core dTrinchieri & Goudfrooiil2002l) . There is also 
a wedge to the south west that is fainter than might otherwise be 
expected. Since no radio emission is associated with the X-ray cav- 
ities, here we conservatively classify them as 'possible' rather than 
'definite'. Deeper, low frequency radio observations are required. 



12 Faint Images of the Radio Sk y at Twenty Centimetres, 

http://sundog.stsci.edu/ dBecker et aljl995h 13 Giant Metrewave Radio Telescope 
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6 POPULATION STUDY 

6.1 Prevalence of Radio Activity 

Our sample contains the X-ray brightest, early-type (elliptical and 
SO) galaxies in the local universe. The sample is drawn from a 
parent catalogue of optically-selected, optical magnitude-limited 
galaxies that is itself approximately 90 per cent complete. Using 
deep VLA observations we have investigated the radio fluxes and 
morphologies of our galaxies. Only one of the 18 galaxies in the 
sample is not detected at radio wavelengths (NGC499). Thus, 94 
per cent of our sample have some level of radio activity at their 
cores. 

The one galaxy for which no central radio source is detected, 
NGC499, is an X-ray bright lenticular (E/SO) galaxy. The upper 
limits on the radio emission from this galaxy are stringent. Com- 
pared to all other galaxies, NGC 499 is very deficient in its 1 .5 GHz 
radio luminosity. 

Of the 17/18 galaxies with detected radio emission, 10 (55 per 
cent) have extended radio lobes. We note that most of our observa- 
tions are only at one frequency, and future observations, especially 
at low frequencies, may detect extended features in some sources 
which so far only exhibit a point source. In all cases (apart from 
IC 310) where the radio morphology is clearly extended, there are 
features in the X-ray emission which correlate with the radio lobes, 
although in one case they are faint. The converse is not quite true. In 
total there are 12 galaxies (67 per cent) which have disturbed X-ray 
morphologies, which in some cases are clearly bubble-like. These 
12 include all nine which have extended radio emission (exclud- 
ing IC 310). The other three (NGC 533, NGC 5044 and NGC 5846) 
have bubble-like features in the X-ray emission, but show no clear 
link to the high frequency radio morphology. 

For the radio point sources, deeper and lowe r frequency ob- 
serva tions with the Long Wavelength Arra y (LWA.lEllingson et ail 
120091) or Low Frequency Array (LOFAR. Ide Vos et alj|2009l ; both 
LOFAR and the LWA are currently under construction.) would be 
useful in searching for past episodes of activity from the central 
AGN. 

Unfortunately 4/18 of our galaxies are observed at 4.9 GH0 
We therefore note that the observations of these galaxies are there- 
fore slightly biased against detecting extended radio emission. 
Hence, the number of galaxies which have extended radio emission 
could be a lower limit on the true number. 

Only two of the eighteen galaxies in our sample, IC 310 and 
NGC 4203, host strong X-ray emitting AGN. However, the inclu- 
sion of these two galaxies does not affect our results or conclusions 
significantly. 



6.2 The link between cool cores and radio activity 

Recent observations of the action of AGN at the centres of groups 
and clusters have alluded to a fe edback process between the AGN 

Jchurazov et al.ll2002l ; lAllen et all 



Table 4. BCES fits to T cool vs L R 



and t he surrounding gas (see e.g. 



l2006t rMcNamara & Nulsen 2007). The cooling of the X -ray emit- 



ting gas provides the fuel for the AGN, whose output is there- 
fore regulated by the amount of gas cooling. However, the output 



14 The VLA archive was searched to find the radio observations which 
showed the most detail for these galaxies, with the aim to find 1.4 GHz 
observations. 4.9 GHz observations were only used when no suitable 
1.4 GHz observations were found. 



Method 


Slope 


Intercept 


Bootstrap Slope 


BCES(Y— X) 


-0.138 ±0.0515 


3.96 ± 1.93 


-0.146 ±0.0631 


BCES(X — Y) 


-0.383 ± 0.147 


13.3 ± 5.56 


-0.457 ± 0.926 


BCES Bisector 


-0.257 ±0.0714 


8.47 ±2.66 


-0.275 ±0.123 


BCES Orthogonal 


-0.142 ±0.0538 


4.14 ±2.07 


-0.151 ±0.0664 



The slopes from the BCES estimator are also shown in Fig. [6] 



of the AGN prevents the gas from cooling as rapidly as it oth- 
erwise would. Estimates of the energy injected by the AGN into 
the central regions of the clusters and galaxies appears to match 
that required to prevent (excessive) cooling (e.g. iBirzan et alj2004l : 
iDunn & Fablanll2006r . iRaffertv et alj fcooe). However, the details 
of the coupling of the injected energy to the intracluster medium 
or galactic halo are not fully understood. Attempts at modelling 
this self-re gulated feedback have made some progress ove r recent 
years (e.g. Vernaleo & Revnolds 2006; Briiggen et al. 2009), but as 
yet there is no observational evidence of the coupling mechanisms. 
Hence, this picture of feedback is conceptually attractive, although 
a number of details remain unclear. 

To test this scenario, we calculate the mean mass-weighted 
cooling times within 1 kpc for all the galaxies and compare them 
to Lr, the monochromatic radio luminosity (see Fig. [6). Here, 
Lr — vbJjv-^ where i/r, and L VR are the observing frequency and 
monochromatic radio luminosity as shown in Table [5] The cooling 
times indicate how rapidly the gas is cooling onto the SMBH, and 
hence give a rough measure of the fuelling rate of the AGN. The 
radio luminosity, as we include the large scale emission from the 
radio lobes, gives a crude measure of the mechanical output of the 
AGN from the relati vistic jets. For mor e precise calculations with 
a smaller sample, see lAllen et alj J2006T) . 

In many of the X-ray observations of the galaxies in the sam- 
ple a central point source was identified and removed. This creates 
a minimum radius for which our annular deprojection can obtain 
values for the cooling time. For these systems, to estimate the cool- 
ing time we fit simple power laws to the density and cooling time 
profiles and extrapolate these inwards to obtain appropriate values 
for the missing inner annular region. Monte Carlo simulations are 
used to estimate the final uncertainties in the results. 

Setting aside the two galaxies for which the X-ray emission is 
dominated by nuclear sources (IC 310 and NGC 4203) we observe 
(see Fig. |6) a slight anti-correlation between the radio luminosity 
and th e mass-weighted coo l ing ti me, with a slope of —0.14. The 
BCES [Akritas & Bershadvl dl996t) slope estimators are shown in 
Table|4] and take into account uncertainties in both parameters and 
the possible presence of intrinsic scatter. Although the sample is 
small and the scatter is large, we do see indication of a link be- 
tween the central cooling times and the radio luminosity which is 

intriguing. 

Our findings are consistent with the results of lAllen et ail 

(2006) and lend further support to the idea that the X-ray gas fuels 
the central AGN in these systems. 

Even in this comparatively small sample, the link between ex- 
tended radio emission and features in the X-ray gas is not totally 
clear. Although all sources which have extended radio emission 
have features in the X-ray gas, the converse is not true. In study- 
ing the effect of AGN on their surroundings, using only the X-ray 
morphology could give misleading results. The selection involved 
in finding AGN-inflated cavities are not presently well understood 
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Figure 6. The mass-weighted cooling time within 1 kpc against the L„ R 
for all galaxies except IC310 and NGC4203. The best fit lines shown are 
the BCES Bisector (solid) along with the errors on the Bisector (dotted) are 
shown. The colour scale shows the mass-weighted central entropy within 
1 kpc for the galaxies. 

dNulsen et ai]|2007l ; iBirzan et alj|2009l) . The fraction of galaxies 
showing AGN infla ted cavities in this sample is around twice what 
iNulsen et alj j2007l) report. The high level of core radio detections 
in the sample presented here indicate that, (a) the fraction of AGN 
currently injecting mechanical energy into their surrounding could 
be much higher than previously thought, and (b) the criteria for an 
AGN to be "active" need to be quantified. 

The AGN in t he elliptical galaxies with detected cavities in 
Nulsen et ail d2007l) appear to be injecting more energy than re- 
quired to offset the cooling of the X-ray gas. This may result if the 
cavities for an average elliptical are very difficult to detect, and so 
only exceptional (overpowered) outbursts are (usually) observed. 
For lower mass systems (groups and ellipticals), the overall level of 
(mechanical) energy injection into their surrounding environments 
remains incompletely understood. 

6.3 Radio-to-X-ray flux ratios 

We have determined integrated radio-to-X-ray flux ratios for the 
galaxies. The measured quantity, GTi-x, follows the radio-to-X-ray 
loudness parameter for nucl ear AGN emission, 7Zx, determined by 
iTerashi ma & Wilson ( 2003). We write 

G1I x = u £i L v JLk. (1) 

where Lx is the 0.1 — 2.4 keV band X-ray luminosity. 
ITerashima & Wilson! I2003T) note that the boundary between radio- 
loud and radio-quiet nuclear sources occurs at log 10 (7?,x) = —4.5. 

The distribution of log 10 (C/7?.x) is shown in Fig. [7] with the 
data detailed in Table \5\ As can be seen from Fig. [7] there is a 
large spread in \og(G1Zx) ■ Also shown are the results of fitting a 
Gaussian model to the distribution: the peak of the Gaussian lies at 
log 10 (Sftx) = -4.02 (-4.08) with a variance of 0.78 (0.79) for 
a downhill simplex (least squares) fit. 

We have also calculated GTZx, nvss values ba sed on the ra- 
dio fl ux densities determined from the NVSS Survey dCondon et al.l 
1998). We show the monochromatic radio luminosities from the 
NVSS, Lr, nvss, along with the C?"ft x , nvss in Tabled The NVSS 
flux densities were all measured at 1.4 GHz using a more compact 




-7 -6 -5 -4 -3 -2 

log 10 (CR x ) 



Figure 7. The distribution of log(Q1Zx ) for the sample of 17 galaxies 
with detected radio emission. The hashed bars are the group dominant- , 
the cross-hashed are the cluster dominant- and the solid bars are the non- 
dominant galaxies. The two fitted curves are for a least-squares (dashed) 
and a downhill simplex (solid) Gaussian model fit to the total distribution. 
IC 3 10 and NGC 4203 are included in this plot. 

configuration of the VLA. More faint emission is included at the 
price of the lower spatial resolution. There are some differences 
between the values of QTZx calculated from the high spatial reso- 
lution VLA observations and those from the NVSS. However, all 
but one (IC 310) differ by less than a factor of two. All galaxies in 
our sample have NVSS detections. See SectionfTJfor more on the 
NVSS properties. 

Our sample contains four dominant cluster galaxies 
(NGC 708, NGC 1399, NGC 4486, NGC 4696), seven domi- 
nant group galaxies (IC 1860, NGC 507, NGC 533, NGC 1550, 
NGC 4406, NGC 4649, NGC 5044) and seven galaxies (IC310, 
NGC 499, NGC 1404, NGC 4203, NGC 4472, NGC 4636, 
NGC 5846) that are not the dominant galaxy of a group or cluster. 
We see no significant difference in the G7lx values for dominant 
and non-dominant galaxies. 



6.4 Prevalence of Radio-Loud AGN 

There is a perception that the fraction of AGN in the general popu- 
lation of galaxies is small, even though SMBH are present in most 
moderately sized galaxies. Part of this perception is h istorical: early 
studies only found AGN fra ctions of a few per cent jDressler et al.l 
ll985l : lHuchra & Burdll992h . However, as the data have improved , 
this fraction has increased dCarter etalj|200ll : iMiller et al.l [20031 : 
S antr a et al .120071) with nuclear activity being found in galaxies that 
would not, classically, have been classed as AGN. For example, in 
the Palomar sampl e of nearby bright gala xies (magnitude limited), 
iNagar et all d2005l) and lFilho etai] d2006l) used high resolution ra- 
dio imaging to show that at least a quarter of galaxies, and probably 
more, are AGN. However, many of these AGN are low luminos- 
ity. In another study, a combined optical and X-ray study of Ellip- 
tical galaxies in the Virgo Cluster find X-ray signatures of AGN 
in the cores of 49-87 per cent target galaxies with stellar masses 
> 10 10 M© dGallo et alj|2008l) . At lower galaxy masses the frac- 
tion falls to 3-44 per cent. This decline in the AGN fraction with 
de creasing galaxy m ass has been seen at the cluster and group level 
bv lBest et al]d2005t) . 
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Table 5. Radio Loudness 



Source Type Frequency L x S„r Lr GTZx £r, NVSS G"Kx,NVSS 







(GHz) 


( er g e 


- 1 ) 


( er g E 




l Hz 


- 1 ) 


( er 


? E 


- 1 ) 










( erg s 


- 1 ) 










IC 1860 


P 


1.365 


5.89 


X 


10 42 


1.39 


X 


10" 


25 


2.12 


X 


10 39 


3.60 


X 


10" 


-5 


2.86 


X 


10 38 


4.86 


X 


io- 


-5 


NGC499 


P 


1.365 


2.24 


x 


10 42 












- 






- 














- 






NGC 507 


E 


1.525 


5.75 


X 


10 42 


8.74 


X 


io- 


25 


8.27 


X 


10 3S 


1.44 


X 


10" 


-4 


5.36 


X 


10 38 


9.31 


X 


io- 


-5 


NGC 533 


P 


4.885 


1.95 


X 


10 42 


8.91 


X 


io- 


26 


2.42 


X 


10 38 


1.24 


X 


io- 


-4 


2.23 


X 


10 38 


1.14 


X 


io- 


-4 


NGC 708 


E 


1.365 


1.23 


X 


10 43 


8.25 


X 


io- 


25 


4.72 


X 


10 38 


3.83 


X 


iir 


-5 


3.85 


X 


10 38 


3.13 


X 


io- 


-5 


NGC 1399 


E 


1.465 


4.90 


X 


10 41 


4.63 


X 


10" 


24 


3.06 


X 


10 38 


6.25 


X 


10- 


-4 


1.31 


X 


10 38 


2.68 


X 


io- 


-4 


NGC 1404 


p 


1.365 


1.78 


X 


10 41 
























2.46 


X 


10 36 


1.38 


x 


10" 


-5 


NGC 1550 


E 


1.365 


7.24 


X 


10 42 


1.27 


X 


io- 


25 


5.07 


X 


10 37 


7.73 


X 


io- 


-6 


7.51 


X 


10 37 


1.04 


X 


io- 


-5 


NGC 4203 


P 


4.885 


1.74 


X 


10 41 


1.80 


X 


io- 


25 


3.18 


X 


10 37 


1.83 


X 


10- 


-4 


3.09 


X 


10 36 


1.78 


X 


io- 


-5 


NGC 4406 


P 


4.885 


1.29 


X 


10 42 


4.91 


X 


io- 


27 


8.35 


X 


10 35 


6.48 


X 


io- 


-7 
















NGC 4472 


E 


1.489 


3.09 


X 


10 41 


2.23 


X 


io- 


24 


1.16 


X 


10 38 


3.74 


X 


io- 


-4 


1.07 


X 


10 38 


3.47 


X 


io- 


-4 


NGC 4486 


E 


1.435 


1.02 


X 


10 43 


8.73 


X 


io- 


22 


7.01 


X 


10 40 


6.85 


X 


io- 


-3 


6.75 


X 


10 40 


6.60 


X 


io- 


-3 


NGC 4636 


E 


1.425 


4.47 


X 


10 41 


6.45 


X 


io- 


25 


3.20 


X 


10 37 


7.16 


X 


io- 


-5 


3.79 


X 


10 37 


8.49 


X 


io- 


-5 


NGC 4649 


E 


1.465 


2.19 


X 


10 41 


2.82 


X 


io- 


25 


1.44 


X 


10 37 


6.58 


X 


10- 


-5 


1.42 


X 


10 37 


6.48 


X 


io- 


-5 


NGC 4696 


E 


1.565 


1.95 


X 


10 43 


3.64 


X 


io- 


23 


1.07 


X 


10 40 


5.50 


X 


io- 


-4 
















NGC 5044 


P 


1.465 


6.31 


X 


10 42 


2.99 


X 


io- 


25 


5.49 


X 


10 37 


8.70 


X 


io- 


-6 


6.09 


X 


10 37 


9.65 


X 


io- 


-6 


NGC 5846 


P 


1.465 


5.13 


X 


10 41 


1.02 


X 


io- 


25 


1.08 


X 


10 37 


2.10 


X 


io- 


-5 


2.12 


X 


10 37 


4.13 


X 


10" 


-5 


IC310 


E 


1.465 


3.98 


X 


10 42 


1.50 


X 


10" 


24 


1.21 


X 


10 39 


3.05 


X 


io- 


-4 


1.30 


X 


10 39 


3.26 


X 


io- 


-4 


NGC 4203 


P 


1.465 


1.74 


X 


10 41 


7.86 


X 


io- 


26 


4.15 


X 


10 36 


3.10 


X 


io- 


-4 


3.09 


X 


10 36 


1.78 


X 


10" 


-5 



The Types show whether the radio morphology is P-point like or E-extended. The Lr are calculated at the frequency shown in the table, without any scaling 
to a common frequency. We also show the £r,nvss calculated from the NVSS, which are all at the same frequency. We separate IC310 and NGC 4203 
which appear to be a different class of source as the other galaxies. 



In determining the fraction of radio l o ud ga laxies present 
in a sample of SDSS galaxies iBest et al] J2005t) use a limit 
for radio loudness of Li.4GHz > 10 23 WHz~ 1 . We converted 
this limit to a luminosity limit at the frequency of the NVSS 
catalogue (1.5 GHz). Therefore, any galaxy with L r, nvss > 
1.4 X 10 39 erg s _1 would be classed as radio-loud in IBest et al.l 
(2005). In our sample only NGC 4486 is above this limit; how- 
ever NGC 4696 would also be if it was in the NVSS catalogue, 
and IC 310 falls directly on the boundary. This gives a radio loud 
fraction of 11 per cent, (2/18) and 17 per cent (3/18) if including 
IC 310. We show the distribution of radio luminosity with X-ray lu- 
minosity in Fig.[8] along with the threshold value for radio-loudness 
from lBest et aLl]2005h . 

This radio loud fraction is much lower than the radio de- 
tected fraction in our sample (17/18). Although using a cut-off in 
the Li.4GHz ensures selecting AGN radio activity out to large dis- 
tances, the results from this study shows that it does not include 
a large number of AGN with low radio luminosities. Of the AGN 
which are not included when using a radio luminosity cut-off, a sig- 
nificant fraction have extended radio emission, though others are 
just radio cores. To obtain a clearer picture of the level of radio 
activity within a sample of galaxies, the radio-to-X-ray flux ratio 
(Section [6.3t of a galaxy can be used. The two galaxies which are 
classed as radio loud from the radio luminosity limit both have ex- 
tended radio lobes. However a further 8 (including IC310) also 
have extended radio emission, indicating active mechanical energy 
injection into their surroundings. Therefore the radio luminosity 
limit results in an incomplete picture of the level of energy injection 
into the surroundings of elliptical galaxies. 

Although our sample selects massive galaxies which are X- 
ray bright, we are not purposefully selecting AGN host galaxies; 
central point sources typically account for only a very small frac- 
tion of the total X-ray flux. Our results indicate that AGN activity 
is common in most massive galaxies; massive galaxies being the 




42.0 42.5 43.0 

loglO-^X, 0.1-2.4 keV 



43.5 



Figure 8. The radio luminosity (Lr,„) as a function of the X-ray luminosity 
(^x.n.l — 2.4kev) along w ith a dashed line showing radio loudness cut-off 
used by Best et al. 12005). The red squares are extended radio sources and 
the blue circles are radio point sources. The purple triangle is NGC 5846 
which although being a point source has clear indications of interactions in 
the X-ray emission. The cyan diamond is NGC 1404 using NVSS data. The 
upper limit for NGC 499 is also shown. The limitations of using a cut-off in 
the Z/i.4 GHz are clear. 



ones which are more likely to have and be able to retain a halo of 

X-ray gas. 

In a large sample of clusters and groups of galaxies, IBest et all 

d2007l) find that brightest group and cluster galaxies are more likely 
to host a radio loud AGN than a field galaxy of an equivalent mass. 
This increased likelihood also extends to galaxies which are close 
to the cluster centre. Not all of our targets are dominant galaxies. 
Most are in clusters or groups. The two most radio loud sources 
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in our study are the brightest galaxies of the Virgo and Centaurus 
Clusters. 



6.5 Particle Content 



Here w e follow the analy sis in iDunn & FabiarJ J2004h:l Dunn et al. 
(2005) and more rec ently [Croston et al. j2005h : |Pe Yound d2006); 
Shurk in et al.l d2008l) and lCroston et alj ( [20080 in studying the par- 
ticle content of the radio lobes. Under the assumption that the rela- 
tivistic radio-emitting plasma is in pressure balance with the X-ray 
gas, we are able to determine the energy, and hence particle content 
of the radio plasma. 

Measurements of the synchrotron emission from the radio 
lobes can be used to calculate the energy contained within the 
relativistic electrons present in the lobes. iFabian et alj (2002) 
studied the lobes of 3 C 84 in the Perseus cluster. Subsequently 
IDunn & Fabianl d2004l) : IDunn et alj d2005h investigated a larger 
number of radio sources in gala xy clusters. For a mo re detailed 
description of the method see e.g. IDunn & Fabianl d2004) . 

For a continuous synchrotron spectrum with a single spectral 
index, a, between v\ = 10 MHz and v-2 = 10 GHz, the energy in 
relativistic elections is 



E c 



a + 0.5 



-B 



-3/2 



erg 



aB 



-3/2 



(2) 
(3) 



Also taking into account the energy within the magnetic field, the 
total energy in the lobes is 



E tot = kE c + Vf^-, 
on 



(4) 



where k accounts for any other particles present in the lobe which 
are not accounted for by the simplistic model spectrum. V = 
^ttRiR^ is the volume of the bubble, where Ri and i? w are the 
bubble radii along and perpendicular to the jet axis respectively. / 
is the volume filling fraction of the relativistic plasma. 

The magnetic field present inside the bubble is estimated by 
comparing the synchrotron cooling time of the plasma to the age 
of the lobe. The latter can be estimated from the sound-speed ex- 
pansion timescalj^l No strong shocks are seen in the X-ray gas, so 
sound speed expansion timescale gives a lower limit on the age of 
the bubbles, and hence an upper limit on the magnetic field. 

We obtain upper limits on the ratio k/ /, 



k 

- < 
f 



Sir J a 



(5) 



for all galaxies with extended radio emission, where the 3V arises 
from the energy density of the relativistic particles. If k/f — 
1, i.e. log 10 (fc//) = 0, then the lobe is filled with a purely 
electron-positron plasma with emission only in the range 10 MHz 
to 10 GHz. Under the assumption / ~ 1, then if k/f > 1, k > 1, 
which implies that there are "extra" particles required for the lobe 
to be in pressure balance with its surroundings. These could be ther- 
mal protons, mixed into the relativistic plasma as the jet travels 
out from the AGN, or they could be electrons which radiate out 
of the assumed region. Since we assume a simple power-law slope 
for the spectral index, any large deviation at low (and hence unob- 
served) radio frequencies will also change the calculated k/f. We 



15 Since the radio lobes of interest are still attached to the active radio jets, 
this is the most sensible, easy to calculate, timescale. 



,£3 

s 

Z 




3 4 
log 10 (k/f) 

Figure 9. The distribution of the upper limits on k/f for the 16 galax- 
ies which have extended radio emission. We show both a standard binning 
(hashed) and a Monte-Carlo routine which was used to account for the large 
uncertainties of the k/f values (solid). The errorbars are from the Monte- 
Carlo routine and are for the solid bars. 



have also determined the value of k/f if the sources were to be in 
equipartition. However in none of the radio bubbles does equipar- 
tition between the relativistic particles and the magnetic field lead 
to pressure balance between the radio and X-ray plasmas. 

As can be seen from Table[6] our sample exhibits a large range 
of k/f upper limits. The distribution of values is also shown in Fig. 
[9] (We use a Monte Carlo algorithm to estimate the uncertainties 
in the upper limits). The dis tribution of k/ f i s broad, in agreemen t 
with the conclusions of e.g. lDunn et al. ( 2005); Bi rzan et"al]d2007h . 



6.5.1 Low Frequency Radio Spectrum 

The value of k/f is expected to increase as a plasma ages with 
synchrotron losses, causing the average energy of the particles to 
drop. Electrons with low energies will emit at frequencies below 
the observed energy range and act as "extra particles" in the model. 
Since bubbles will only appear 'young' for a short time, most will 
be seen when k/f is large. Eventually the radio emission will drop 
below detection limits, leading to an upper limit to the values of 
k/f estimated from GHz radio data. There is a suggestion from Fig. 
[9] that the distribution peaks at around 10 3 , but with asymmetric 
tails either side. 

We have investigated whether, by modifying the spectrum to 
include additional low frequency radio emission from low energy 
particles (extending v\ = v m i n < 10 MHz), it is possible to 
achieve k/f equal to one. The minimum possible radio frequency 
is the cyclotron frequency as estimated from the magnetic field es- 
timates. For around half of the bubbles k/f — 1 is achieved before 
the frequency limit is reached. The bubbles for which k/f — 1 is 
not achievable are those which have the largest k/f values to start 
with. 

In some of the galaxies in our sample, there is no extended 
1.4 GHz radio emission, even when there are features in the X- 
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Table 6. k/f Results 



Galaxy 


Lobe R\ 


H W 




S 


a 


log(B) 


log(fc//oqptn) log(fc// sound ) log(^ min ) 


log(fc//„ mi J 




(kpc) 


(kpc) 


(kpc) 


(Jy) 











NGC507 


E 


10.70 


15.20 


14.90 


0.027 


-1.40 ±0 


20 


-5.09 ±0.02 


2 


83 ± 0.34 


2.77 ±0 


34 


1 


04 ± 0.95 


-0.00 ± 0.02 


NGC507 


W 


11.00 


14.30 


12.20 


0.048 


-1.40 ±0 


20 


-5.10 ± 0.02 


2 


54 ± 0.35 


2.47 ±0 


34 


4 


37 ±0.87 


-0.00 ±0.02 


NGC708 


E 


4.93 


2.72 


5.78 


0.014 


-0.60 ±0 


20 


-4.86 ± 0.03 


.3 


52 ± 0.17 


3.34 ±0 


16 


1 


58 ± 0.03 


2.68 ± 0.78 


NGC708 


W 


4.75 


2.82 


4.75 


0.017 


-0.60 ±0 


20 


-4.85 ± 0.03 


3 


44 ±0.18 


3.28 ±0 


IS 


1 


59 ± 0.03 


2.61 ± 0.78 


NGC1399 


N 


2.88 


1.94 


8.39 


0.11 


-0.90 ±0 


10 


-4.73 ± 0.03 


2 


56 ± 0.17 


2.62 ±0 


17 


1 


67 ±0.03 


0.01 ± 0.02 


NGC1399 


S 


3.96 


2.50 


8.42 


0.13 


-0.90 ±0 


10 


-4.81 ± 0.03 


2 


60 ±0.18 


2.65 ±0 


17 


1 


63 ± 0.02 


0.01 ± 0.01 


NGC1550 


E 


4.34 


6.37 


5.60 


0.0035 


-1.00 ±0 


50 


-4.83 ± 0.03 


6 


83 ± 0.58 


6.63 ±0 


57 


1 


62 ± 0.03 


-0.01 ±0.02 


NGC1550 


W 


7.34 


4.55 


10.70 


0.0029 


-1.00 ±0 


50 


-4.97 ±0.03 


6 


51 ± 0.57 


6.25 ±0 


57 


1 


46 ± 0.03 


-0.01 ±0.02 


NGC4472 


E 


2.24 


2.24 


4.81 


0.029 


-0.80 ±0 


10 


-4.68 ± 0.03 


3 


65 ± 0.15 


3.70 ± 


15 


1 


72 ± 0.03 


2.08 ± 0.60 


NGC4472 


W 


3.07 


1.83 


2.82 


0.038 


-0.80 ±0 


10 


-4.77 ±0.03 


3 


27 ±0.14 


3.30 ±0 


11 


1 


66 ± 0.02 


1.66 ± 0.61 


NGC4486 


J 


1.46 


0.91 


1.49 


35.8 


-0.60 ±0 


20 


-4.47 ±0.02 


1 


13 ± 0.22 


0.97 ±0 


22 


1 


96 ± 0.02 


0.01 ± 0.02 


NGC4486 


CJ 


1.29 


1.29 


1.48 


56.5 


-0.60 ±0 


20 


-4.44 ±0.02 


1 


23 ± 0.16 


1.11 ±0 


16 


1 


99 ± 0.03 


0.01 ± 0.01 


NGC4636 


NE 


0.71 


0.41 


0.82 


0.028 


-0.60 ±0 


20 


-4.40 ±0.02 


2 


43 ± 0.18 


2.59 ±0 


18 


2 


04 ± 0.02 


0.01 ± 0.02 


NGC4636 


SW 


0.77 


0.39 


0.86 


0.0185 


-0.60 ±0 


20 


-4.42 ± 0.01 


2 


61 ± 0.16 


2.76 ± 


16 


2 


02 ± 0.01 


0.01 ± 0.02 


NGC4649 


N 


0.51 


0.49 


0.84 


0.0029 


-0.60 ±0 


20 


-4.25 ± 0.02 


4 


13 ± 0.16 


4.26 ±0 


16 


2 


17 ±0.01 


3.64 ±0.70 


NGC4649 


S 


0.65 


0.42 


1.03 


0.0019 


-0.70 ±0 


20 


-4.32 ± 0.02 


.3 


93 ± 0.22 


4.06 ± 


21 


2 


12 ± 0.02 


3.04 ±0.95 


NGC4696 


E 


2.83 


1.75 


2.32 


1.3 


-0.75 ± 


20 


-4.72 ± 0.02 


1 


53 ± 0.25 


1.40 ± 


25 


1 


68 ± 0.03 


0.01 ± 0.01 


NGC4696 


W 


2.88 


1.96 


3.24 


0.65 


-0.75 ± 


20 


-4.73 ±0.02 


1 


94 ±0.24 


1.81 ±0 


21 


1 


68 ± 0.03 


0.01 ±0.01 



R\ and R„ are the dimensions of the bubbles with Rdist their separation from the AGN. fc// eq n is the equipartition value of k/f and fc// soun< j 
is the value as estimated from the sound speed expansion timescale. u m i n is the lower cut off to the radio spectrum required for k/f = 1 if it can 
be reached before the cyclotron frequency. k/f VlDin is the value of k/f at that frequency. The bubb le descriptors are N-Nort h, S-South, E-East, W - 
West, J- Jet CJ-Counter Jet. The spectral i ndices have been taken fr om the follo wing: NGC 507 IColla et al J [l975h. NGC708 faardcastle et alj Eoblh . 
NGC 1399 iEkersetalJll989h. NGC4649 IStanger & Warwickl [l98r3). NGC4486 iHines et alj l l989h . NGC 463 6 IStanger & ^ Warwick! 1 1986k ~NGC 4472 
(High -resolution X-ray and radio observations of the giant elliptical galaxies NGC 4636 and 4649IVolh . NGC 4696 iFabian et alj2005h " 



ray emission which point towards an interaction (e.g. NGC 5044 
and NGC 5846). Studies which have included low fr equency ra- 
dio emission (e.g. lDunn et alj2005l : lBirzan et alj2 008) obtain k/f 
values which are larger than the values from GHz radio emission. 
However, the increase in k/f for ghost bubbles, regardless of the 
type of radio emission is currently detected from them, is at least in 
part an effect of not being able to sample the full radio spectrum of 
the lobes. 

As the electrons age, the synchrotron radio flux in the GHz 
band falls. Therefore more "extra particles" are required for pres- 
sure balance, resulting in a higher value calculated for k/f (as- 
suming / = 1). However, the aged synchrotron electrons are still 
present in the lobe and still supporting it, but radiating at much 
lower frequencies. Therefore a measure of k/f at GHz frequencies 
does not fully describe the k/f of the complete lobe. Observations 
with the new low frequency ra dio telescopes (LOFAR, LWA or 
the Murchison Widefield Array Lonsdale & MWA Collaboration 
l2007t) will uncover the spectral shape at low frequencies and so 
allow further investigations into the detected and the presumed cav- 
ities. 



6.5.2 Possible Biases 

In the calculations outlined above, we have used the sound speed 
timescale to estimate the age of the bubbles, and compared this 
age to the synchrotron cooling time of the plasma to estimate the 
magnetic field. These are not the only possibilities available when 
calculating the k/f for the cavities. 

Older cavities have been described by a buoyancy timescale, 
the time taken to rise buoyantly in the potential well of the galaxy 
( Churaz ov et al. 2000). This buoyant rise takes place after the initial 
expansion, and separates the cavity from the galactic nucleus. In all 



of the cavities studied here, they are still attached to the nucleus, 
with on the whole, radio emission also connected to the radio core. 
The buoyancy timescale is therefore unlikely to be a reasonable 
description for the age of the cavities. 

Another timescale of interest is the t ime for the intra-galacti c 
medium to refill the displaced volume jMcNamara et al. 2000). 
However, this is also most applicable to older cavities which have 
detached from the nucleus and risen up, with the IGM flowing 
around behind the cavity. 

No strong shocks are seen around the cavities at the centres of 
nearby clusters of galaxies (e.g. Perseus, A2052), the expansion of 
these bubbles is not currently highly supersonic. Although fed by a 
relativistic jet, and hence at some point expand supersonically, on 
average the sound speed is an appropriate value to calculate the age 
of the cavities from. We assume that this is also the case in these, 
in some cases lower-power, cavity systems. 

To obtain an estimate on the magnetic field of the relativistic 
plasma, we compare the synchrotron cooling time to the age of the 
cavity. Using the observation that GHz radio emission is detected, 
we conclude that the synchrotron cooling time must be longer than 
the age of the cavity. However, we are unable to say how much 
longer; we can only place a lower bound on the synchrotron cooling 
time, and hence an upper bound on the magnetic field strength. 
Future low-frequency radio observations will help in determining 
the spectral shape of the radio emission, and hence improve the 
estimates on the magnetic field strength, and hence the age of the 
bubbles. 

Improving the age estimates of the bubbles may be possible 
with increasingly accurate simulations of the behaviour for the evo- 
lution of radio lobes in clusters of galaxies. For older bubbles, a 
combination of the sound speed expansion and the buoyancy rise 
time may be the best current description for the ages. 
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7 EXTENDED SAMPLE 

We have attempted to extend our study to lower X-ray and radio 
fluxes. Dropping the flux cut to 1 x 10~ 12 erg s _1 cm" 2 increases 
the sample from 18 galaxies to 42. Of these 42 galaxies, 25 (60 
per cent) hav e their L x values deter mined from follow-up pointed 
observations JO'Sullivan et alj|2001f) . In order to limit our study to 
the most massive elliptical and SO galaxies, we introduce a lumi- 
nosity cut along with the flux cut. Limiting the galaxies to those 
with Lx > 10 erg s _1 ensures we do not select very nearby 
low mass galaxies. 

To determine the radio properties of these galaxies we use the 
NVSS survey at 1.4 GHz for the northern hemisp here, and the Syd- 
ney University Mo longlo Sky Survey (SUMSS. iBock et alj[l999l ; 
iMauch et alj2003l) at 843 MHz for the southern hemisphere. These 
surveys, albeit at slightly different frequencies, are comparable in 
terms of their sensitivity and resolution. Using these, we are able 
to determine whether the galaxy hosts a radio source. We perform 
a search using a radius of 30 arcsec in NVSS and then a 60 arc- 
sec search for those galaxies where nothing is found. We use the 
SUMSS postage stamp cutouts to check those galaxies which are 
out of the NVSS survey area. 

Table [7j summarises the radio information from the NVSS 
and SUMSS for all galaxies in the extended sample. It is clear 
that this extended sample still has a very high radio detection rate, 
with 34/42 (81 per cent) galaxies having a radio detection within 
NVSS/SUMSS. The radio sensitivity of this analysis is much lower 
than that in the previous sections, using pointed VLA observations. 
Therefore the total fraction could be significantly higher. In Fig.llQI 
we show the distribution of galaxies and those which have radio 
detections in the NVSS and SUMSS. 

For the extended sample we are not able to comment on the 
morphology of the X-ray emission, and only very coarsely on the 
radio morphology. However, we can use the classifications given to 
the galaxies in the SIMBAD3 database, which are noted in Table 
[7] There are no galaxies with classifications which would indicate 
that they have radio emission that we have overlooked/missed with 
the two surveys used. The galaxies which do not have any radio 
emission are a mixture of field galaxies, as well as those in groups 
or clusters, but not dominating them. 

We investigated whether from the SIMBAD morphological 
markers we could clearly determine whether the X-ray emis- 
sion is likely to be dominated by an X-ray bright AGN. How- 
ever NGC4472, NGC4486 and NGC4636 have markers of either 
Seyfert 2 galaxies or a LINER type AGN but all have clear and 
strong diffuse X-ray emission in addition to the AGN emission. 
Without high spatial resolution X-ray observations of these galax- 
ies we cannot comment on the ratio of the nuclear to the diffuse 
X-ray emission in these galaxies. 



8 LUMINOSITY FUNCTION 

Fig. [TT] shows the cumulative radio luminosity functions for the 
main and extended galaxy samples. The data > 1 x 10 1 have been 
fitted with a simple Schechter function, 



Table 7. Extended Sample 



axy 


ux 
m y 


urvey 


Type 


K. CUM 111 I 


^X 

'°glO er S S 


IC310* 


168.1 


N 


GiC 


0.0156 


42.60 


IC1860* 


18.3 


N 


GiC 


0.0221 


42.77 


IC2006 


- 


N 


GiC 


0.0045 


<41.09 a 


IC4296 


546.6 


N 


rG 


0.0118 


41.59 


IC4765 


32.4 


S 


GiG 


0.0144 


41.89 


NGC57 


_ 


N 


G 


0.0136 


41.71 


NGC410 


5.8 


N 


GiG 


0.0140 


41.97 


NGC499* 


- 


N 


GiG 


0.0136 


42.35 


NGC507* 


61.7 


N 


IG 


0.0166 


42.76 


NGC533* 


28.6 


N 


GiG 


0.0157 


42.29 


NGC708* 


65.7 


N 


rG 


0.0136 


43.09 


NGC741 


478.8 


N 


rG 


0.0151 


41.79 


NGC777 


7.0 


N 


Sy2 


0.0136 


42.14 


NGC1399* 


208.0 


N 


GiC 


0.0045 


41.69 


NGC1404* 


3.9 


N 


GiC 


0.0045 


41.25 


NGC1407 


87.7 


N 


GiG 


0.0051 


41.06 


NGC1550* 


16.6 


N 


G 


0.0120 


42.86 


NGC2300 


2.9 


N 


IG 


0.0069 


41.22 


NGC2305 


- 


S 


G 


0.0114 


41.73 


NGC2329 


363.7 


N 


rG 


0.0175 


42.18 


NGC2340 


_ 


N 


G 


0.0182 


42.14 


NGC3091 


2.5 


N 


GiG 


0.0126 


41.69 


NGC4073 


_ 


N 


GiG 


0.0195 


42.44 


NGC4203* 


6.1 


N 


LIN 


0.0040 


41.24 


NGC4261 


4066.7 


N 


LIN 


0.0078 


41.27 


NGC4406* 


_ 


N 


GiC 


0.0040 


42.11 


NGC4472* 


219.9 


N 


Sy2 


0.0040 


41.49 


NGC4486* 


138487.0 


N 


LIN 


0.0040 


43.01 


NGC4636* 


77.8 


N 


LIN 


0.0040 


41.65 


NGC4649* 


29.1 


N 


GiP 


0.0040 


41.34 


NGC4696* 


5674.0 


S 


GiC 


0.0092 


43.29 


NGC4936 


39.8 


N 


LIN 


0.0102 


41.75 


NGC5044* 


34.7 


N 


GiG 


0.0075 


42.80 


IN uLj UyU 


1 jZj.o 


C 
.3 


lu 


U.UlUJ 


41. JJ 


NGC5129 


7.2 


N 


GiG 


0.0224 


42.20 


NGC5328 




N 


GiC 


0.0152 


41.94 


NGC5419 


349.2 


N 


rG 


0.0132 


41.84 


NGC5846* 


21.0 


N 


GiP 


0.0057 


41.71 


NGC6407 


54.6 


S 


G 


0.0145 


42.00 


NGC6868 


138.8 


S 


GiG 


0.0086 


41.29 


NGC7049 


93.2 


s 


GiC 


0.0068 


41.07 


NGC7619 


20.3 


N 


GiG 


0.0099 


41.69 


Total 


34 


42 
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exp 
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* - those galaxies in the detailed sample. The type-codes are as follows 
G-Galaxy, GiC-Galaxy in Cluster, GiG-Galaxy in Group, GiP-Galaxy in 
Pair, rG-Radio Galaxy, IG-Interacting Galaxy(ies), Sy2-Seyfert 2 galaxy, 
LIN-LINER-type Active Galaxy Nucleus. The surveys used are N-NVSS 
and S-SUMSS. a IC 2006 only has an upper limit on the X-ray luminosity 
so whether it is truly part of this sample is debatable. It also has no radio 
detection, and so it is possible that, given the high radio detection rate, that 
it has a much smaller X-ray luminosity than shown in the table. 



A best fit results in f = 0.0289, a = -0.367 and L, = 7.79 x 
10 24 W/Hz for the bins > 1 x 10 21 . Below around ~ 2 x 10 21 W/Hz 
the function breaks again and can be modelled with a powerlaw of 

slope —0.1. 

We also use the model used bv lBestetal.l ( l2005h 
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Figure 10. The di stance-luminosity distribution of the galaxy sample of 
Beui ng et alj 19991 with the lower flux cut and the luminosity cut used to cre- 
ate the extended sample. Those with detected radio emission in the NVSS 
or SUMSS have a red circle, those with no detection have a crossed,empty 
circle. We show the distance and X-ray luminosity limits, as well as the old 
(dotted) and new (solid) flux limits used to create this sample. 
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L yvss ,W/Hz 

Figure 11. The luminosity functions of the main (red squares) and extended 
(blue circles) samples using the NVSS radio luminosities, supplemented by 
the SUMSS catalogue in the extended sample. We also show the best fit 
relation to the extended sample using both a Schec hter function (gre en line) 
and the form of the luminosity function shown in Best et al. 1 2005) scaled 
for a mass of 1O 8 ' 66 M0. as the solid black line. 

where / = 3.5 x 10" 3 , a = 1.6, /3 = 0.37, 7 = 1.79 and 
L* = 3.2 x 10 24 W/Hz but scaled for a mass of 1O 8 66 M . In our 
study we have not split the sources by black hole mass, as in this 
case both samples are too small for this to be appropriate. Unsur- 
prisingly the luminosity function is not very smooth, especially for 
the main sample presented in this paper. Despite this, the behaviour 
of the luminosity function of the radio galaxies shown in lBest et al.1 



(2005), with the break at around 10 24 W/Hz, fits the luminosity 
function well. 

When using both models, there appears to be evidence for a 
second break around 2 x 10 21 W/Hz, and the behaviour below this 
break appears to show that the fraction continues to rise towards 
one. This secon d break, if true, o ccurs at lower luminosities than 
those present in lBest et alj (2005), but they do see a flattening of 
the luminosity function for the case of the most massive black hole 
bins. 

Given the close agreement bet ween the shapes o f the luminos- 
ity functions presented here and in iBest et al.l |2005), the galaxies 
in our samples do not appear to be special - at least in their radio 
luminosities. Furthermore, given the comparatively few galaxies in 
our sample (42) compared to the many more in the SDSS study 
(2215) and the good agreement between the two shapes of the radio 
luminosity function, this implies that our sample is fairly unbiased 
to the radio luminosity of the galaxies within it. 



9 CONCLUSIONS 

Using a parent catalogue of elliptical galaxies with well determined 
optical properties we have identified a sample of 18 nearby, X-ray 
bright elliptical galaxies. Out of these 18 galaxies, 17/18 (94 per 
cent) have detected central radio emission and over half (10/18) 
have extended radio structures larger than 1 arcsecond. 12 galaxies 
have clear evidence for some disturbance of their X-ray halos, nine 
of which also have extended radio emission indicating a likely in- 
teraction of the radio and X-ray plasmas. Only two of t he galaxies 
in our sample are above the radio loudness criterion of IBest et al.1 
(2005) (2/18, 1 1 per cent). However, the true level of radio activity 
is much higher. Extending the sample to lower X-ray fluxes, at the 
expense of data quality and completeness we find that 34/42 galax- 
ies have a radio detection with the NVSS and SUMSS surveys. The 
lower sensitivity of these two surveys means this percentage should 
be regarded as a lower limit. 

The mass-weighted cooling time within 1 kpc appears to anti- 
correlate with the total radio luminosity, particularly once systems 
for which the X-ray emission is dominated by nuclear point sources 
(IC 310 and NGC 4203) are excluded. This lends further support to 
the idea that hot gaseous halos of elliptical galaxies provide an im- 
portant fuel source for the central AGN. 

We have calculated upper limits on k/ f, where k is the ratio 
of the total particle energy to that of relativistic electrons radiating 
in the range of 10 MHz and 10 GHz and / is the volume filling 
factor of the plasma in the cavity, for the extended radio lobes in 
the sample. The distribution of the upper limits on k/f is broad and 
consistent with what would be expected as a radio plasma ages. 

We calculate an X-ray radio loudness parameter for the total 
galactic emission in these two bands, QTZx- There is no clear cor- 
relation of GTZx with the galaxy environment or Lx- 
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